Conservation of energy in coherent backscattering of light 
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Although conservation of energy is fundamental in physics, its principles seem to be violated in the 
field of wave propagation in turbid media by the energy enhancement of the coherent backscattering 
cone. In this letter we present experimental data which show that the energy enhancement of 
the cone is balanced by an energy cutback at all scattering angles. Moreover, we give a complete 
theoretical description, which is in good agreement with these data. The additional terms needed 
to enforce energy conservation in this description result from an interference effect between incident 
and multiply scattered waves, which is reminiscent of the optical theorem in single scattering. 
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Conservation of energy is one of the most fundamental 
principles in physics. Not only in mechanics, but also 
in thermodynamics it rests on fundamental symmetries, 
which so far have not been violated. There have however 
been instances, where novel effects gave rise to suspicions 
of non-conservation of energy as in the distribution of 
positron energies in beta decay. This was only explained 
by the introduction of the neutrino in Fermi's four-point 
theory [1]. Another instance which seemingly contradicts 
the conservation of energy is the coherent backscattering 
cone which appears when waves propagate in turbid me- 
dia. In that case, a twofold enhancement of the backscat- 
tcrcd intensity with respect to Lambert's law for diffuse 
scattering is observed. This enhancement decays over an 
angular scale of (kl*) , where k is the wavenumber of 
the wave and I* is the transport mean free path of the 
medium [2]. Due to its fundamental nature, it can be ob- 
served in areas as different as solid state physics [3], soft 
matter physics [4] , astro- [5] and geophysics [6] , and with 
various kinds of waves like sound waves [7] , microwaves 
[8] and, as in the experiments presented in this letter, 
visible light [4, 9]. 

The additional energy contribution of the coherent 
backscattering cone to the intensity can however not be 
explained by a corresponding reduction in transmission 
at surfaces not considered in the experiment. Rather, the 
backscattering cone is also observed from samples which 
can theoretically as well as experimentally be treated as 
filling an infinite half-space, meaning that the waves can 
reach no other surface than the one considered [2] . More- 
over, different polarization channels also show a backscat- 
tering cone [10], such that the energy in the cone cannot 
be obtained from another polarization channel. 

The origin of the backscattering cone lies in the inter- 
ference of waves propagating along reciprocal paths. This 
interference can only spatially re-distribute the backscat- 
tered energy. Thus the energy enhancement at small an- 
gles has to be accompanied by a corresponding energy 
cutback to ensure conservation of energy. However, such 
an energy cutback has so far not been observed experi- 
mentally. Moreover, it is not described by the prevailing 
theories [2, 10, 11]. This can be problematic as the scal- 



ing of the width of the backscattering cone with kl* is 
commonly used to characterize multiple scattering ma- 
terials. In particular in turbid samples, when the cone 
becomes very broad, there has to be a sizeable correc- 
tion if conservation of energy is to hold. This is of great 
importance in the study of Anderson localization of light 
[12], where a reliable knowledge of the parameter kl* is 
needed to characterize the phase transition from diffusive 
transport to an insulating state [13, 14]. 

In this letter, we present measurements of coherent 
backscattering, where the incoherent background is de- 
termined on an absolute scale. With this we are able 
to show that there is a reduction in backscattering in- 
tensity at all angles compensating for the enhancement 
in the back-direction. A theoretical description which 
fits these data shows that the reduction in enhancement 
results from a new interference effect between the inci- 
dent and multiply scattered waves. This is analogous to 
the shadow term, which accounts for flux conservation 
in the optical theorem [15]. In multiple scattering, the 
terms needed to ensure energy conservation correspond 
to the so-called Hikami box or quantum crossing [16, 17]. 
The present experiment constitutes a direct observation 
of a scattering process described by a Hikami box, which 
plays a central role in quantum mesoscopic physics [17]. 

Our main setup to study the angular distribution of 
the backscattered light consists of 256 photosensitive 
diodes attached to a semi-circular arc with a diameter 
of 1.2 m, in the center of which the sample is located 
[18]. In this way, we can detect light over a range of 
-60° < < 85° with a resolution of 0.14° for |0| < 10°, 
~ 1° for 10° < \6\ < 60° and ~ 3° for 9 > 60°. For the il- 
lumination a continuous wave dye laser with a wavelength 
of 590 nm is used. The measurements are done using cir- 
cularly polarized light in order to reduce the influence of 
single scattering. To average over random speckle pat- 
terns, the samples are rotated. As the very tip of the cone 
at 8 ~ can not be resolved with this setup, the central 
part of the backscattering cone, \6\ < 3°, is measured 
separately using a beam splitter and a charge-coupled 
device (CCD) camera to a resolution of 0.01° [18]. 

In order to calibrate the photodiodes, we have used a 
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FIG. 1: Backscattering cone of R700 evaluated ignoring 
(dashed curve) and taking into account (full red curve) ab- 
sorption in the reference sample. The former is positive for 
all angles, resulting in an uncompensated energy enhancement 
of the cone. For the latter, energy enhancement and energy 
cutback are balanced. This is quantified by the half-space in- 
tegral of the intensity / = —0.005(7). The inset shows a time 
of flight measurement of the teflon reference, which follows 
diffusion theory with an absorption time of 3.3 ns. 

block of teflon as a reference sample. Teflon has a trans- 
port mean free path of ~ 300 /im and hence the backscat- 
tering cone of teflon at a wavelength of 590 nm has a 
FWHM of about 0.02°, which is much narrower than the 
angular resolution of the wide-angle setup. This implies 
that with the wide-angle setup teflon can be considered 
as a purely incoherent signal, which is properly normal- 
ized given by fi(j + fi/(fi + 1)) [17], where fi = cos (9 and 
7?* describes the distance over which the diffuse inten- 
sity enters the sample. These diode signals are measured 
at several different incident laser powers, which are de- 
termined independently with a calibrated power-meter. 
Interpolation of the measured data then yields a calibra- 
tion function for each photodiode [18]. To be able to 
determine the intensity of the backscattered cone on an 
absolute scale, the incoherent background needs to be 
known. Since the cone of teflon cannot be resolved in the 
wide-angle setup, the teflon reference measurements can 
be considered to describe the incoherent background of 
the TiC>2 samples. In doing so however, one neglects the 
different albedo of teflon with respect to the sample. The 
proper background for the sample, ai llc , samp is therefore 
given by that of the reference, ai nCjro f multiplied by the 
ratio of the albedos of sample and reference, A samp /A re {. 

For an estimation of the albedos to a level of better 
than one percent, one needs to take into account losses 
at the sample/reference boundaries, as well as losses due 
to absorption. Up to now these losses have not been fully 
taken into account in the evaluation of the backscattering 
cones [18, 19]. 



FIG. 2: Measurements of backscattering cones for different 
samples. As the cone width increases, more energy needs to 
be compensated. Thus the most turbid samples (R700 and 
Aldrich) lead to a noticeable energy cutback at angles around 
45° . The amplitude of this cutback is reduced with decreasing 
turbidity, but stays positioned around 45°. This indicates 
that the cutback is due to effects occurring at a fixed length 
scale close to A. 

The loss due to leakage is estimated by comparing the 
diffuse energy in an infinite half-space at a certain time 
with the amount of energy that is left in a volume-cutout 
of the infinite half-space corresponding to the size of the 
reference [11]. This will neglect edge effects, but since 
we are only dealing with the very tails of the diffusive 
cloud, this will be accurate to the desired level. For this 
estimate, the diffusion coefficient of the sample/reference 
needs to be known. The loss factor for absorption can be 
calculated from the integral over the path length distri- 
bution P(D, t, t), where D is the diffusion coefficient and 
r is the absorption time. 

Both quantities needed to calculate the albedo, D and 
r can be determined with a time resolved transmission 
experiment [14, 20], see the inset of Fig. 1. Such a time 
of flight experiment directly gives the path length distri- 
bution inside the sample, which is a function only of D 
and r. In our experiment, the same dye laser is used as 
in the backscattering experiments, thus making a possi- 
ble wavelength dependence of D and r irrelevant. We 
thus calculate the ratio of the albedos using the absorp- 
tion length L a = \J Dt with respect to the reference and 
hence the absolute level of the incoherent background for 
all samples. Subtracting this background then directly 
gives a proper measure of the backscatter enhancement. 

To represent cones with a wide variety of cone widths, 
we have used samples of ground TiC>2 particles in its ru- 
tile structure with different particle diameters (R700: 245 
nm and Aldrich: 540 nm), a mixture of TiC>2 (R700) and 
ground chalk in a weight ratio of one to five, as well as 
solid teflon. The TiC>2 particles are commercially avail- 
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able as pigment for white paint [14] . 

As can be seen in Figs. 1 and 2, a cutback of the 
backscattered energy is indeed observed when taking into 
account the different loss factors of reference and sam- 
ple. This is most significant for samples with very wide 
cones like R700 and Aldrich with FWHM w 3.8° and 3.4° 
respectively, where the enhancement is noticeably below 
zero for a range of 50° on both sides of the cone, as shown 
in Fig. 2. The Ti0 2 -chalk-mixturc with FWHM w 1.5° 
still shows a slight lowering of the enhancement, while 
the enhancement for teflon with FWHM w 0.013° is es- 
sentially zero away from the cone. Note that unlike the 
coherent backscattering peak, the energy cutback is not 
characterized by a specific angle but it is rather spread 
out over the whole angular range. Furthermore, these en- 
ergy cutbacks do compensate the energy enhancements 
of the cones, as the integral of the enhancement a{ff) over 
the backscattering half-space I = J a(8) sin 6 d6 is zero 
for all investigated samples within the margins of error. 
This shows that a determination of the absolute intensity 
scale is crucial for the correct observation of the backscat- 
tering cone. For such a determination, the different loss 
factors of reference and sample had to be accounted for. 

Theoretically, the coherent backscattering cone has 
been described in great detail [17]. In the geometry of a 
semi- infinite medium of section S (see Fig. 3a), this de- 
scription makes use of the following well-known expres- 
sion for the coherent albedo a A , 

"° = 4^ / dr i dr 2ff A (r l! r 2 )P(r 1 ,r 2 ) (1) 

where H A (ri,r 2 ) = V 2 e *(k+k')-(ri-r 2 ) f or an m _ 

cident plane wave normal to the interface. P(ri , r 2 ) is the 
probability of having a multiple scattering path starting 
at ri and ending at r 2 . The first factor in H A describes 
the attenuation of an incident plane wave over a distance 
of the order of the elastic mean free path I*. The second 
factor in H A accounts for interference and leads to the 
enhancement with an angular width of order l/kl*. 

The interference term H A is the product of four am- 
plitudes describing the two incoming and the two out- 
going plane waves. It is known as a quantum cross- 
ing and it is at the origin of coherent effects in quan- 
tum mcsoscopic physics such as weak localization, uni- 
versal conductance fluctuations and eventually it leads 
to the localization transition. Energy (or number of par- 
ticles) conservation imposes constraints on the quantum 
crossings. It is well-known [16] that in order to ful- 
fill this constraint, two other contributions H B ' C (ri,r 2 ) 
must be added to H A , which mix the in- and outgo- 
ing wave vectors. Energy conservation thus imposes that 
/ dR (H A + H B + H c ) = 0, where R = n — r 2 . 

A complete description of coherent backscattering 
must therefore include these additional contributions H B 
and H c , which are equal. The physical basis of these 
contributions lies in a coupling of the light fields at the 
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FIG. 3: Wave configurations of the contributions H A and 
H B C to coherent backscattering. H A (a) describes interfer- 
ence between time reversed amplitudes (full and dashed lines) 
and gives the classical cone shape [2]. When ri and V2 are 
within a transverse distance of A, the time-reversed loops have 
to be considered as closed and the amplitudes are becoming 
coupled. This is described by H B ' C (b), as an interference 
effect between the incident plane wave and the attenuated 
spherical wave traveling between the points ri^ and a newly 
introduced scatterer r at an angle (3 (see text). 

first and the last scatterer (ri and r 2 ), when they are 
within a volume of order A 2 1*. This coupling originates 
in an interference of the incoming plane wave with the 
multiply scattered spherical wave and is described by in- 
troducing an additional scattering event located in r (see 
Fig. 3b). This is reflected by the short range behavior of 
i? B (ri,r 2 ) and the additional contribution results from 
almost closed diffusive trajectories. Consequently, this 
interference is not restricted to small angles 9 as is the 
case for a A . For the case of an incident wave normal to 
the interface, the contribution a B can be written as: 

a c~^J dvP{v,v)e~^h\ (2) 

where 

/ ikr' ■]* 

dt , e ,..,_ e -r,^_ (3) 

is calculated to leading order in (fcZ*) -1 . It is interest- 
ing to note the similarity between h and the shadow 
term which occurs in the optical theorem and ensures 
flux conservation for single elastic scattering. Here, the 
shadow terms h describe the interference between the in- 
cident and the multiply scattered waves. The integral 
in (2) can be solved approximately to give a correction 
a B+c _ -i.i5(kl*)- 2 (j,/(n + 1). Thus the correction 
is of order — (kl*)~ 2 . Noting that the angular integral of 
a A ~ (kl*)~ 2 , we indeed retrieve the energy conservation 
condition of quantum crossing, namely that the integral 
over a A+B+c is zero, as it should be. 

The expression for a B+c shows that the interference 
effect in (2) is twofold. First, h is purely imaginary so 
that h 2 is negative resulting in a depletion of the coherent 
albedo proportional to (kl*)~ 2 . Secondly, this interfer- 
ence term does not contribute at a specific angular value 
but it is rather spread out over the whole angular range. 
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FIG. 4: Comparison of the backscattering cone of R700 with 
corrected and uncorrected enhancement. The agreement be- 
tween the measured data and the fit of the corrected enhance- 
ment a^ +B+c (dashed) is perfect within the errors. The un- 
corrected enhancement (dotted) calculated with the value 
of kl* obtained by the fit of the data with a^ +B+c describes 
the cone itself quite well, but shows significant deviations in 
the area of the energy cutback. 



ference between incoming and multiply scattered waves 
at short distances. Since the cone is basically the Fourier 
transform of this intensity this leads to a broadly dis- 
tributed energy cutback [10]. This improved description 
of the cone-shape for extremely turbid samples also al- 
lows a reliable determination of kl* in such samples. A 
comparison of these result with those obtained previously 
on the same samples show however that the values of kl* 
thus determined change only very little. This implies 
that the systematic dependence of deviations from clas- 
sical diffusion is indeed as it was previously determined 
[13]. 
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A fit of the total coherent albedo af +B+c to our data 
is shown in Fig. 4 by the dashed line. The dotted line 
shows only a^, where the same value of kl* has been 
used as in af +B+c . Comparing this description with 
the uncorrected data in Fig. 1 shows that the kl* values 
determined with those data are close to their real values. 

In order to properly determine the turbidity of the 
sample, the value of kl* obtained from (1) still has to 
be corrected for internal reflections [21], which lead to 
a broader distribution of light paths, and thus to a nar- 
rower cone. Furthermore, the ratio between the transport 
and elastic mean free paths will influence the pre-factor 
in the correction but is assumed to be close to unity. 

In conclusion, we have shown experimentally that co- 
herent backscattering docs fulfill conservation of energy. 
To this purpose, the losses of the reference sample had to 
be quantified via a time of flight measurement to ensure 
an absolute energy calibration of the setup. If the loss 
of the reference differs significantly from the loss of the 
sample, this leads to different positions of the incoherent 
background in spite of equal incident laser energies. 

Furthermore, we have provided a complete theoretical 
description of coherent backscattering based on the cal- 
culation of the three terms H ' B ' C that contribute to the 
Hikami box or quantum crossing. H describes the steep 
angular variation around backscattering and H B,C take 
into account crossed diagrams dressed by a scattering 
impurity. Such an impurity provides an additional inter- 
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